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Abstract. We present the bundle (Aff(3) ® C <g> A)(IR 3 ), with a geometric 
Dirac equation on it, as a three-dimensional geometric interpretation of the 
SM fermions. Each (C <S> A)(M 3 ) describes an electroweak doublet. The Dirac 
equation has a doubler-free staggered spatial discretization on the lattice space 
(Aff (3) (g> C)(Z 3 ). This space allows a simple physical interpretation as a phase 
space of a lattice of cells in IE 3 . 

We find the SM SU(3) C x SU{2) L x U(1) Y action on (Aff (3) ®C ® A)(M 3 ) 
to be a maximal anomaly-free special gauge action preserving E(3) symmetry 
and symplectic structure, which can be constructed using two simple types of 
gauge-like lattice fields: Wilson gauge fields and correction terms for lattice 
deformations. 

The lattice fermion fields we propose to quantize as low energy states of a 
canonical quantum theory with Z2-degenerated vacuum state. We construct 
anticommuting fermion operators for the resulting Z2-valued (spin) field the- 
ory. 

A metric theory of gravity compatible with this model is presented too. 

1. Introduction 

After the success of relativity, the interest of modern physics has been centered 
on four-dimensional spacetime. If a concept requires a preferred frame, this is, for 
many physicists, already sufficient to reject it. This seems unjustified, given that in 
the most popular approach — string theory — general-relativistic symmetry is only 
an effective, derived symmetry, and the fundamental Minkowski space structure is 
as unobservable as a preferred frame. 

Of course, a theory with preferred frame has to explain the observable relativistic 
symmetry. But this is possible: In [25], GR in harmonic gauge - and, especially, 
the Einstein equivalence principle - is derived from axioms for a condensed matter 
theory. The basic ideas of this derivation we present here in appendix [X] Once 
this problem is solved, there seems to be no decisive argument against a preferred 
frame. 

One of the assumptions of this theory of gravity is that matter fields are 
not exterior fields, but describe material properties of the same medium. As a 
consequence, to obtain a complete theory, we need a condensed matter model 
for the SM fields too. The aim of this paper is to present such a model. 
We have found a three-dimensional geometric interpretation of SM fermions as 
(Aff(3) ® C <£> A)(]R 3 ), together with a doubler-free staggered discretization on the 
lattice space (C <£> Aff(3))(Z 3 ), which allows a condensed matter interpretation as 
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the phase space of a lattice of cells located in R 3 . Moreover, this model allows, 
essentially, to compute the SM gauge group and its action on the fermions. Thus, 
all SM fields observed until now can be described in this way. 

Let's start with the bundle (Aff(3) <g) C <g> A)(M 3 ), which we propose as a three- 
dimensional geometric interpretation of the SM fermions. The bundle (C Cg) A)(M 3 ) 
describes an electroweak doublet J*] Each of the 3 • (3+ 1) components (a^) 6 Aff(3) 
of an affine transformation we associate with such an electroweak doublet: The 
upper index i denotes the generation, /i = the leptonic sector, \x > the quark 
sector, where the three positive values n S {1, 2, 3} define the three quark colors. 

On the bundle (C (B> A)(]R 3 ) exists a three-dimensional geometric Dirac operator 
- an analogon of the Dirac-Kahler operator [TI5] on (C (£> A)(]R 4 ). This operator is 
sufficient to define the Dirac matrices a 1 . We find also natural operators Ii as well 
as f3 = 7° on (C ® A)(IR 3 ). The Dirac equation we define in its original Dirac form 
idtip = Hip, as an evolution equation on (C (8 A)(M 3 ). This equation contains eight 
complex fields and describes a doublet of Dirac particles. 

In analogy with the staggered discretization of the four-dimensional bundle 
(C (£> A)(R 4 ), we have also a staggered discretization of the three-dimensional Dirac 
operator. It lives on a three-dimensional spatial lattice 1? . It is a staggered dis- 
cretization, with only one complex component on each lattice node, and eight differ- 
ent types of lattice nodes. Similar to the four-dimensional staggered discretization 
of A(M 4 ) on Z 4 (see [21 El SI El [24l [26] ) , it is a doubler-free discretization of the Dirac 
equation on (C ® A)(R 3 ). In other words, we obtain a lattice evolution equation on 
a three-dimensional lattice C(Z 3 ), which gives, in the continuous limit, two Dirac 
fermions. 

For all SM fermions (the bundle (Aff(3) <g) C <g) A)(M 3 )) we obtain a first order 
lattice equation on (C ® Aff(3))(Z 3 ). This lattice space allows a physical interpre- 
tation as the phase space for a three-dimensional lattice of elementary cells, where 
the state of each cell is described by a single affine transformation from a standard 
reference cell (see figure [TJ. 

This physical interpretation gives us two important structures: First, a symplec- 
tic structure of the phase space, second, a natural action of the Euclidean group 
E(3). These structures may be used to restrict the action of the gauge group. 
For a compact gauge group, we can always construct a preserved Euclidean metric 
on the phase space, which, together with a preserved symplectic structure, allows 
to construct a preserved complex structure. Thus, preservation of the symplectic 
structure requires unitarity of the gauge groups. 

The left action of E(3) on Aff(3)(Z 3 ) transforms the lattice as a whole. The 
requirement of preserving this symmetry for the gauge groups consists of two parts: 

• To commute with rotations, gauge groups have to preserve generations and 
to act on all three generations in the same way. This holds for all SM gauge 
fields. 

• To commute with translations, one direction in the leptonic sector has to 
be preserved. All SM gauge fields leave right-handed neutrinos and their 
antiparticles invariant, thus, a common invariant direction exists in the SM. 



Independent of this paper, three-dimensional geometric fermions have been proposed by 
Daviau 9 , and the idea that geometric fermions may be used to describe electroweak doublets 
has been proposed by Hestenes I16| . 
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Figure 1. The space (C ® Aff(3))(Z 3 ) of the lattice model sug- 
gests an interpretation as the phase space (with configuration space 
Aff(3)(Z 3 )) of a three-dimensional lattice of deformable three- 
dimensional cells. The configuration of each cell is described by 
an affinc transformation from a standard reference cell. 



Thus, for an appropriate identification of the invariant direction, all SM gauge 
fields preserve E(3) symmetry. 

The lattice theory also leads to another important restriction for the gauge fields: 
We have to define an appropriate lattice model for the gauge fields. A well-known 
way to put gauge fields on the lattice are Wilson gauge fields. Their modification 
to a three-dimensional lattice with continuous time is straightforward. But Wilson 
gauge fields cannot act in a nontrivial way inside the doublets (C <g> A)(M 3 ), because 
these are represented on the lattice as C(Z 3 ), which leaves U(i) as the maximal 
possible Wilson gauge field. Thus, Wilson gauge fields have the same charge on all 
parts of a doublet. The maximal group of Wilson gauge fields compatible with this 
restriction, E(3) symmetry, and symplectic structure is U(3) = SU(3)qcd x U(1)b- 

We find another modification of the lattice equations which, in the large distance 
limit, leads to a gauge-like interaction term for fermions. It describes correction 
terms for lattice deformations^] The coefficients depend only on the geometry 
of the lattice, thus, the resulting equation preserves doublets (C <E> A)(]R 3 ), and 
acts identically on all doublets. We construct generators as operators commuting 
with the lattice Dirac equation, which correspond to the basic lattice shifts on 
the staggered lattice. By themself, they give the gauge group U(2)l x U{2)r. 
One choice of the maximal group of this type, compatible with E(3) symmetry, is 
generated by chiral U(2) l = SU(2) l x U(1)l, and a vector field U(1)q with charge 



The idea that gauge fields may be obtained as effective fields in condensed matter theory to 
describe various types of lattice defects is quite old, see, for example, 1 181 ■ 
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The EM field does not fit into any of the two types. But it can be constructed as 
a combination of them, by the simple formula Q — 21b + Q- Thus, our two types 
of gauge-like lattice fields are already sufficient to construct all SM gauge fields. 

Based on arguments from condensed matter theory, we argue that the gauge 
group should be special, that means, that the charges have zero trace. This excludes 
the field U(\)b, which would, otherwise, be allowed. What remains as the maximal 
possible gauge group G max is the SM gauge group Gsm, together with only one 
additional gauge field — an "upper axial gauge field" , with charge Ijj = j 5 (I 3 — 
which acts on upper quarks and leptons. This extension of the SM is anomalous, 
thus, if we add, as a last conditions, anomaly freedom, we obtain Gsm as a maximal 
anomaly-free subgroup of G max . Thus, we have, essentially, computed the SM gauge 
group from first principles. 

Of course, there are a lot of things left to future research. We have not considered 
yet the Higgs sector and the mass terms. They break the left E(3) symmetry 
preserved by the gauge action. Thus, to describe them, we need some spontaneous 
_E(3)-symmetry breaking. Once the broken symmetry is E(3), it is not clear if we 
need a separate Higgs sector at all. This has to be left to future research. 

What about quantization? The first problem is fermion quantization. We use 
classical, commuting c-number fields in the lattice theory (C (£> Aff(3))(Z 3 ), not 
Grassmann variables as in the Berezin approach to fermion quantization. Thus, we 
need a completely different quantization scheme for fermions. 

We propose such an alternative in section [5] We consider canonical quantiza- 
tion of a field with Z2-symmetric degenerated vacuum state. The lowest energy 
states of such a field define a Z2-valucd (spin) field, yet with commuting opera- 
tors on different lattice nodes. Fortunately, the algebra of (anticommuting) lattice 
fermion operators and the algebra of (commuting) spin field operators appear to 
be isomorph. Unfortunately, the isomorphism is nonlocal and depends on an order 
between the lattice nodes. But we can motivate a particular choice of this order. 
For this choice, we obtain the staggered lattice Dirac operator (exactly in dimen- 
sion one, approximately in higher dimensions) from a much simpler, non-staggered, 
direction-independent spin-field Hamilton operator, which we have obtained from 
canonical quantization. 

Our lattice gauge fields appear as effective fields. Their quantum effects are 
similar to phonons in condensed matter theory. Therefore, they are "quantized" 
already by the quantization of the fundamental (high energy) degrees of freedom 
fields — the states of the elementary cells (and the material between them). In- 
stead of a separate quantization procedure for these fields, the effective quantum 
properties of these fields have to be derived from the fundamental quantum theory. 
The details have to be worked out by future research. Nonetheless, some differences 
are obvious. Especially, different gauge-equivalent configurations describe different 
physical states. Therefore, the standard rejection of anomalous gauge fields, which 
is based on the BRST quantization procedure, which requires strong gauge invari- 
ance, seems premature. That means, the additional field U{\)u may as well appear 
in the final theory. 

The gravitational field is, similarly, an effective large distance field, which de- 
scribes density, average velocity, and the stress tensor of the medium. As well, 
for these fields there will be no separate quantization procedure, but, instead, a 
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derivation of quantum effects of these effective fields from the fundamental quan- 
tum theory. This has been left to future research too. 

2. Geometric interpretation of SM fermions 

Let's consider at first the geometric interpretation of the SM fermions as sections 
of the bundle (Aff(3) <g> C <g> A)(E 3 ). Throughout this paper, we do not consider 
mass terms. Without the mass terms, the three generations of SM fermions appear 
completely identical, simply as three identical copies of the same representation of 
the SM gauge group SU(3) C x SU(2) L x {7(1). 

The group Aff(3) is the group of three-dimensional affine transformations y l = 
djxi + a l on ]R 3 . Each et^ we can identify with an electroweak doublet of the SM 
according to the following simple rules: The upper index i, 1 < i < 3 defines the 
generation. The translational components g,q we identify with the leptonic sector. 
The linear part aj, j > we identify with the quark sector. The lower index j, 
1 < j < 3 denotes the color of the quark doublet. 

This identification of the 3 x (3 + 1) SM doublets with a 3 x (3 + 1) affine 
matrix may be considered, up to now, as pure numerology. But it defines a natural 
action of the Euclidean group E(3), by multiplication from the left. This action 
commutes, as we will see, with all SM gauge fields and plays an important part in 
our computation of the SM gauge action. 

Each electroweak doublet is defined by the bundle (C <8> A)(]R 3 ). ft is assumed 
here that right-handed neutrinos exists, so that neutrinos form usual Dirac parti- 
cles. Thus, qualitatively there is no difference between electroweak quark doublets 
and electroweak lepton doublets: Above contain two Dirac particles. The bun- 
dle (C (£> A)(M 3 ) consists of three-dimensional complex inhomogencous differential 
forms 

K,e{o,i} 
= ^ooo (x) 

^ + i>ioo(x) dx 1 + ipoi (x) dx 2 + V>ooi(a:) dx 3 

+ ipiio(x) dx 1 A dx 2 + ipon(x) dx 2 A dx 3 + V^ioi^) dx 1 A dx 3 
+ ipm(x) dx 1 A dx 2 A dx 3 . 

Thus, we have 1 + 3 + 3 + 1 = 8 complex functions, which gives two Dirac 
fermions. This allows a physical interpretation in terms of a standard model elec- 
troweak doublet. The use of a three-dimensional bundle is essential: In Minkowski 
spacetime, we have only the bundle C x A(M 4 ), with the Dirac-Kahler equation [19] . 
But this equation describes four Dirac fermions. 

On the external bundle A(]R d ) exists a natural geometric Dirac operator D as a 
square root of the Laplace operator A = D 2 . For a general metric, the definition is 
given in appendix [B| In the Euclidean case g^ u — 5^, this Dirac operator has the 
form 

(2) D = d+d* = -ia i d i . 

with operators a 1 , which fulfill the anticommutation relations {a 1 , a J } = 2<5 y . Now, 
together with the skew-symmetric a 1 , it is useful to consider also corresponding 
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symmetric operators ft 1 . They may be defined analogically as 

(3) d-d* = -i(3 i d i . 
Together, they define a set of generators of M 2 d (R) = C/ d ' d (R): 

(4) {a\a j } = 2S ij , {a\p j } = 0, {P\f3 3 } = -25 ij . 

For d = 3, the explicit representation of the matrices a 1 , (3 l is defined by: 



(5) -iafdi^ 
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The last Dirac operator 7 can be obtained now as 

(7) /3 = 7 = fr^fPa^oF-o? = a 1 /? 1 " 2 /? 2 " 3 ^ 3 , 

and appears to be a diagonal operator, which measures the Z 2 -graduation of A(M 3 ). 
The matrices a 1 , (i define a representation of the Dirac algebra 

(8) K,a'} = 2<^; {a i ,0} = O ] (a'f = 1 = 1. 

For the (massless) Dirac equation we prefer to use the original form, as proposed 
by Dirac, with the operators a 1 : 

(9) id t ty = -ia l d^. 
The operators /, defined by 

(10) 2ie llk h = P(3 k , 

define a representation of the isospin algebra su(2). We identify them with the 
(weak) vector isospin Ii = t 1 l + t' 1 r . The Jj commute, as they should, with the Dirac 
equation as well as with 7 . Thus, the operator I 3 may be used to split the bundle 
(C Cg) A)(M 3 ) into two parts with eigenvalues ^3 — ±|, so that each of the parts 
contains a full representation of the Dirac algebra. 

An interesting question is how the spinor representation a lj — a ! a J on the Dirac 
particles is connected with the representation so (3) of geometric rotations of the 
bundle (C <S> A)(R 3 ). The answer is that geometric rotations are generated by the 
operators oj^ defined by 

(11) ^ = tfo? - P 1 ^ = a ij - 2ie ljk I t . 
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Thus, the true, geometric rotations of our geometric interpretation are a combina- 
tion of spinor rotations and isospin rotations. The operator 7 s = — ia 3 a 2 a 1 turns 



out to be the (modified) geometric Hodge * operator ( 90 ) 



2.1. Symplectic structure. We have a complex structure in our geometric in- 
terpretation. Now, every complex structure defines a natural symplectic structure 
uj = dz A dz. We know that all the SM gauge groups are unitary groups, thus, they 
preserve the complex structure. As a consequence, they also preserve the symplectic 
structure. Therefore, we can postulate the following: 

Postulate 1. All gauge fields preserve the symplectic structure derived from the 
complex structure of (C ® A)(]R 3 ). 

The question we want to consider here is if we really need the complex structure. 
May be the symplectic structure is already sufficient? Or do we obtain, in this way, 
some additional gauge fields? No, at least as long as we consider only compact 
gauge groups. For compact gauge groups, we have the invariant Haar measure 
dfi(g), and it has a finite norm. This allows to construct, for a given action of a 
compact group, an invariant Euclidean norm (., .). All we have to do is to start 
with an arbitrary norm (., .) and to compute the average of the Haar measure: 

(12) (a, b) = J (ga,gb} dfi{g). 

The resulting Euclidean distance (., .) is already preserved by the gauge group 
action. Once we have a preserved Euclidean metric (., .) together with a preserved 
symplectic structure u>, we can already construct a preserved complex structure by 
the rule 

(13) w(a,ib) = (a,b). 

As a consequence, our postulate [I] is sufficient to restrict the gauge group to an 
unitary group. Thus, in the geometric interpretation we can forget about the 
complex structure and restrict ourself to the symplectic structure, that means, we 
can interpret the space (Aff(3) <g> C <g> A)(K 3 ) as a phase space. 

2.2. Euclidean symmetry. On Aff(3), we have a well-defined left action of the 
Euclidean symmetry group E(3) C Aff(3). 

The action of the rotation group 0(3) C E(3) extends immediately to 
(Aff(3)®C<g>A)(IR 3 ) as 



(14) ' 



In terms of our interpretation, these rotations rotate the three generations of the 
SM. Now, all SM gauge groups preserve generations. (Remember that we consider 
here the massless case, thus, define generations not in terms of mass eigenstates, 
but in such a way that they contain electroweak doublets completely.) Moreover, 
they act on the different generations in exactly the same way. As a consequence, 
they commute with the action of our group of rotations 0(3). 

Let's extend now the action of the subgroup of translation T 3 C E(3) on 
(Aff(3) (X) C (8) A)(M 3 ). For this purpose, we have to define a shift operator 

(15) t : % -> r(f )%, 
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where r(t) : ^> — > defines a scalar shift operator on (C ® A)(IR 3 ). This shift is 
an additive action of R on (C ® A)(]R 3 ), and it should not depend on x. Therefore, 
it is uniquely defined by a single shift vector c— (c K ) G C as 

(16) c = t(1)*-*, 

which we name the "direction of translation". After this, translations are defined 
as r(t)ip K — > tpK + tc K for all n, and we have extended the definition of translations 
from Aff(3) to (Aff(3) <g> C ® A)(M 3 ). 

In our interpretation, translations act, by shifts, only on the leptonic doublets. 
Once we already have found that rotations commute with all gauge groups, it would 
be nice to have a similar property for translations too. So, what does it mean for 
the gauge groups to commute with translations? The answer is simple — the gauge 
groups have to leave the translational direction c, which is located in the leptonic 
sector, invariant: 

(17) [9,r(t)}=0 «■ gc = c 

Now, the leptonic sector of the SM contains directions which are left invariant by 
all SM gauge fields — the right-handed neutrinos and their antiparticles. Thus, if 
we identify the direction of translation c in such a way, that it is inside the right- 
handed neutrino sector, then all SM gauge fields preserve translational symmetry 
too. 

Thus, with an appropriate definition of the direction of translation c, all SM 
gauge fields preserve the complete E(3) symmetry. This property of the SM gauge 
fields we use in the following as a postulate: 

Postulate 2. All gauge fields preserve the E(3) symmetry defined by the left action 
ofE(3) on .\JjU). 

Note that this observation gives our E(3) symmetry large explanatory power. It 
explains why all SM gauge fields preserve generations and act in the same way on 
the three generations. Moreover, it excludes a lot of very interesting natural and 
symmetric extensions of the SM: 

• The extension of SU(3) C to SU(4) C with lepton charge as a forth color, 
which is part of the Pati-Salam extension of the SM [21] , 

• the left-right-symmetric extension of SU(2)l X U(l)y to U(1)b-l x 
SU(2)l x SU(2)r, which is also part of the Pati-Salam extension of the 

SM EH, 

• and all GUTs which use at least one of these extensions as a subgroup, 
especially 50(10) GUT. 

Indeed, all these extensions of the SM act on right-handed neutrinos in a non- 
trivial way, and, moreover, they also leave no other direction invariant. As a con- 
sequence, they cannot commute with any choice of the direction of translation c, 
thus, cannot commute with Euclidean translations. 

Nonetheless, these principles are not yet sufficient to compute the SM gauge 
group. There remain interesting nontrivial extensions like SU(5) GUT |13j or 
chiral color with SU(3) L x SU(3) R instead of SU(3) C [Hj. 



A CONDENSED MATTER INTERPRETATION OF THE SM 



9 



3. The lattice Dirac operator 

Let's consider now a discretization of our Dirac equation in space, leaving time 
continuous. Using naive central differences, we obtain the following lattice equation: 

(18) idtMn) = Y,- t ( a T^(^'(n+h l )-^(n-h l )) 

i 

on the lattice space f2 = C 8 (Z 3 ). 

It is easy to see that this lattice equation contains eight doublers. Indeed, let's de- 
fined eight so-called "staggered" sublattices, labelled by A = (Ai,A 2 ,A3) € {0, l} 3 , 
defined by the condition 

(19) n x = {ip K (n)\n = k + A mod 2} 

so that f2 = J2>J} X - K is easy to see that the naive lattice Dirac equation (in our 
representation (|5j) preserves the decomposition into the staggered sublattices. As 
a consequence, it is easy to get rid of the doublers, and sufficient to preserve only 
one of the eight sublattices £! 000 , with A = {0,0,0}. Thus, our staggered sublattice 
is defined by the condition 

(20) n = k mod 2. 



This doublcr-frcc lattice equation (18 1, (20 1 can be obtained from a much more 



genereal, geometric construction, which is presented in appendix |B.1| It is the 
same geometric construction, which gives, in the case of the four-dimensional Dirac- 
Kahler equation [19] on the spacetime bundle A(M 4 ), the Kogut-Susskind staggered 
fermions [3j [20j |26] in lattice gauge theory (see [2] 14] 15] Ho] [24] ) . 

Now, it is interesting to see what happens with the other operators we have 
defined in the continuous limit. On f2, the operators a* 7 5 and the shift 
operators r, : ^f(n) — ► ^(n + hi) are well-defined. Unfortunately, they do not 
preserve the decomposition into staggered subspaces. Fortunately, there are natural 
replacements for these operators, which already preserve Q 000 . For the generators 

f3 i of Cf' d (R) we can use: 

(21) a l = a l n, p l = P l r t . 

For the other operators /i,7 5 we can use the same formulas we have used in the 
continuous limit to compute them: 

(22) 7 5 = -icPtfa 1 = 7 5 r 1 r 2 T 3 

(23) 2ie ijk Ii = P 3 j3 k = f3 3 !3 k T 1 T k 

Now, the operators 7 5 and Ii generate an interesting group A of operators associated 
with lattice shifts: 

Theorem 1. The group A of operators generated by 7 5 and 2/j has the following 
properties: 

• It preserves the staggered subspaces f2 A . 

• It preserves the massless lattice Dirac equation. [^] 

• There exists an epimorphism ir : A — * Z 3 named "underlying shift opera- 
tor". 



■^Note one advantage of using the original form D = a z di of the Dirac equation: 7 5 does not 
anticommute, but commute with the (massless) Dirac equation. 
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• Ker7r = {±1} and acts by pointwise multiplication. 

• Va S A : a 2 = (na) 2 . 

For a shift operator r £ Z 3 , the equation 7r(f ) = r defines the operator f modulo 
its sign. 

In the continuous limit, the subgroup of A generated by even shifts becomes 
irrelevant. The corresponding factorgroup consists of the following operators: A — 
{±1, ±2/j, ±7 5 , ±2/ i 7 5 }. This set of operators generates the Lie algebra of the 
group U(2)l x U(2)r, which plays an important role in the following. Especially 
it contains the weak gauge group SU(2)l. 

3.1. The cellular lattice model. Let's forget, for some time, about the staggered 
character of the lattice Dirac equation. Then, the lattice space of the discretization 
O 000 is simply C(Z 3 ), with a single complex number on each lattice node. For all 
SM fermions, we obtain the lattice space (C® Aff(3))(Z 3 ). 

Note also that we have a first order lattice equation on it. This suggests an 
interpretation of (C ® Aff(3))(Z 3 ) as a phase space of some physical system: 

(24) ^(n) = rt(n) + iw;(n). 

with configuration variables f^n) : Z 3 — > Aff(3) and momentum variables 7r^(n). 
On the phase space (C ® Aff(3))(Z 3 ) we have the standard symplectic structure 

(25) W =^^(n)Adi;(n), 

Now, the configuration space Aff(3)(Z 3 ) allows a natural interpretation as a 
regular lattice of deformable cells in R 3 (see figure [lj : The state of each cell is 
described by an affine transformation from a standard reference cell. This reference 
cell is assumed to be located in the origin. 

Now, to have such a simple model is, of course, nice and beautiful. But is it only 
an otherwise useless toy, or is it helpful to explain the physics of the SM? We want 
to show here that this model has physical importance. 

First, of course, this model gives the symplectic structure, which we have used 
in section [2.1 1 to derive unitarity. Thus, the model allows to explain our postulate 

m 

But it seems helpful to explain Euclidean symmetry too. Of course, Euclidean 
symmetry is not a property of the full SM, where the mass matrices break this 
symmetry Thus, we need some spontaneous symmetry breaking to explain the 
SM masses. Nonetheless, the lattice model allows to answer the following simple 
question: Why do we have to use the left action of E(3) on Aff(3), instead of the 
right or adjoint action? Abstract group theory remains silent about this. Instead, 
for a lattice of deformed cells, we can look what happens with the lattice if we apply 
the different actions of E(3): 

We consider an almost regular lattice. Then we have approximately 

(26) ^(n)w<55, 

Now, we see that the left action of a rotation rotates the lattice as a whole, including 
the shifts </?o( n )- Instead, the right action of a rotation leaves the cells on their places 
nih and rotates them around these places. This, obviously, modifies the geometric 
relations between neighbour cells. Only the left action rotates the lattice as a whole, 
leaving the local geometry unchanged. Thus, the left action is much more likely to 
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be a symmetry of the theory — something which cannot be seen otherwise. In this 
sense, our cellular model is useful to explain our postulate [2] as well. 

But the most important consequence of the cellular lattice model is, that we 
can apply now condensed matter theory. Especially we can, in the large distance 
limit, define density, velocity, and a stress tensor, and postulate continuity and 
Euler equations. But this is what we need to incorporate gravity into the model. 
A metric theory of gravity with GR limit, based on such an "ether concept", has 
been proposed in [25) . We give a short introduction in appendix [A| 

4. Lattice gauge fields 

While our postulates [T] and [2] impose strong restrictions for the gauge group of 
the SM, we are yet far away from computing the SM gauge group. There are, yet, 
gauge groups much larger than SU(3) C x SU(2)l x U(1)y, which are compatible 
with our postulates. 

But the consideration of the lattice theory allows to impose another type of 
restrictions: It should be possible to "put the gauge action on the lattice". We will 
see that this gives the additional restrictions, which we need to compute the SM 
gauge group almost exactly. 

4.1. Strong fields as Wilson gauge fields. The classical way to incorporate 
gauge fields into a lattice theory are Wilson gauge fields. The classical formalism 
of Wilson gauge fields, even if it was developed for spacetime lattices Z instead of 
our lattice of cells Z 3 , needs only a sufficiently obvious, minor modification. This is 
caused by the fact that we have no discrete structure in time direction. Formally, it 
looks like time remaining continuous. This requires a mixed form for the definition 
of the gauge field: The temporal component Ao(n,t) € g is, like in the continuous 
case, a function with values in the Lie algebra g, and defined on the lattice nodes. 
Instead, the spatial (vector potential) part Ai is described, as usual for Wilson gauge 
fields, by Lie group valued functions U(n, i,t) G G located on the edges n, n + hi of 
the lattice. The most important, defining property of the Wilson gauge field remains 
unchanged too: The lattice gauge symmetry is defined by a gauge-group- valued 
lattice function g(n) : Z 3 — *■ G, which acts pointwise on the lattice (C <g> Aff(3))(Z 3 ) 
and is uniquely defined by a gauge action G x C x Aff(3) -^Cx Aff(3). The gauge 
transformation acts in the following way: 

(27a) -> (g(n,mi(n,t), 

(27b) A (n,t) -> g(n,t)A (n,t)g~ 1 (n,t) - (9 t .g(n, t))g _1 (n, i), 

(27c) U(n,i,t) -> g(n,t)U(n 7 i,t)g~ 1 (n+h i ,t). 



This definition of the gauge action (27 1 shows, that not all imaginable gauge actions 
may be defined in this way. Indeed, the gauge action can act only on the generation 
and color indices. Inside a doublet (C (8 A)(R 3 ), it can act only in a very restricted 
way: An electroweak doublet with fixed generation i and color fj, is represented on 
the lattice as a lattice field C(Z 3 ), so that there is only a single complex number 
i/jpin) in each lattice node. The only possible Wilson gauge action on the lattice 
C(Z 3 ) is, obviously, an action of U(l), which means, that we obtain the same charge 
on all parts of the doublet. 

Now, this already allows to compute the maximal possible Wilson gauge action, 
which is compatible with our postulates [T] and [2j It should be an unitary group. 
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It acts on all generations in the same way, and preserves the generations, thus, 
does not act on the generation index i. Then, it acts with the same charge on all 
parts of electroweak doublets, thus, cannot act on the doublet indices k. Thus, 
it can act only on the remaining index fi. This gives U{4) as the maximal gauge 
group. Moreover, to commute with translations, it has to leave the translational 
direction c in the leptonic sector invariant. But, because it has the same charge 
on all parts of the leptonic doublets, it has to act trivially on the whole leptonic 
sector [i = 0. What remains is the group C/(3) acting on the color index /i > 0. Its 
special subgroup SU(3) can be, obviously, identified with the color group SU(3) C 
of the SM. What remains is the diagonal U(1)b with the baryon charge Ib- 

Having found an upper bound, let's consider the question if these Wilson gauge 
fields will appear, in some natural way, in our condensed matter model. For this 
purpose, let's assume that there is some other material between the cells. Inhomo- 
gencitics of this material will influence the cells, thus, lead to some modification of 
the equation for the cells. In some approximation, the material located between 
two neighbour cells will influence only those parts of the equation which connect 
these two cells. Then, a Wilson gauge field corresponds to such an influence which 
may be compensated by a modification of the state and momentum of one of the 
neighbour cells. It seems reasonable to expect that such influences of the material 
between the cells appear. 

Thus, the consideration of Wilson lattice gauge fields has given us, almost ex- 
actly, an important part of the SM gauge group the strong interactions. 

4.2. Correction terms for lattice deformations. While the consideration of 
Wilson gauge fields is a sufficiently trivial modification of standard Wilson gauge 
fields, the incorporation of weak interactions requires a non-standard approach 
to lattice gauge theory. This new approach is far away from being completed. 
Nonetheless, it already gives a nice correspondence between the properties of the 
gauge groups which may be, in principle, obtained in this way, and the gauge group 
of the SM. 

Assume our lattice I? is not exactly regular but slightly deformed. This requires 
also a modification of the lattice Dirac equation. What can be said about the 
general form of the corresponding correction terms? 

First, a deformation of the lattice is certainly no reason to use different lattice 
equations for the different components 4']i{ n ) € (C <8> Aff(3))(Z 3 ). However de- 
formed the lattice, the correction coefficients are of geometric nature: They depend 
only on the geometry of the deformed lattice. Thus, the deformed lattice equation 
will be an equation on the same bundle (C ® A)(IR 3 ), independent of the generation 
and color indices i and \i. This leads to the following 

Thesis 1. Correction terms for lattice deformations preserve doublets (C® A)(M 3 ) 
and act on all doublets in the same way. 

But this is a signature of weak forces. Thus, it seems natural to postulate a 
connection between weak interactions and correction terms for lattice deformations. 

Let's consider, therefore, possible correction terms in more detail. We start 
with a regular lattice C(Z 3 ), with a function tp(n) on it, and an undistorted lattice 
Dirac equation. For a slightly deformed lattice, the value in a regular lattice node 
x l = TTijh is no longer ip(m), but has to be interpolated using all neighbour nodes. 
Thus, we correct now every occurrence of tp(m) by a weighted sum over values 
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ij)(m + £) on neighbour nodes (including the node m itself): 

(28) ^{m)->iP{m)+^Al(n)ip{m + t) 

e 

with some set of geometric coefficients Aj-(n). These coefficients depend, in general, 
on the basic node n of the lattice equation containing a term with ^(m), on the 
direction of the neighbour £ £ Z , and, moreover, on the occurrence p of the term 
i/)(m) in the undistorted Dirac equation for node n. Using the lattice shift operator 
t c : t/>(m) — > ip(m + £), we can rewrite the expression as 

(29) (l + £ i«(n)r 5 ) ^(m). 

S 

Now, instead of the lattice shift operators t£, which do not commute with the Dirac 
equation, we prefer to use another set of operators associated with lattice shifts, 
namely the operators f £ e A of theorem [T] which commute with the Dirac equa- 
tion. Fortunately, this is possible, it requires only a redefinition of the coefficients 
A|(n) — * Ap(n): We can replace the t£ by t£ = o^(ti)t£, with o^(n) € C, and put 
the coefficients c^(rt) into the A|(n) with A|(n) = A* (n)(o*(n)) -1 . Using these 
pseudo-shift operators T£ 6 .A, together with the geometric coefficients A|(n), gives 

(30) </>(m) ^ (l + E ^W f e) V'(m)- 

We have seven occurrences p € {0, i±} of "0 in the lattice Dirac equation of node 
n, which gives 

(31) 

idti&Z + Ai(n)(f^An) = ((<#' + ^U»)fc)£>«"(« + K) 

Introducing the denotations 

(32) 4(n)=-i(I? + (n) 4(n) = -<0 t i§(n), 

(33) Af(n) = -i l -{A\ + {n) + if_(n)), Ig(n) = -iig(n), 
this becomes 

z(^' + l i«(n)(f c )«> t ^H = 

(a 1 ):' (-i(<#' + if W(f e )£') (> K "(n + ft,) - iMn - /*)) 

^ +Al(n)(f 6 tf ^" {n + hi)+ 2 ^'' {n - hi) ) 

+ Ai(n)(r e )^ K ,(n) 

Now, there are two straightforward simplifications, which we can use, once we are 
interested only in the large distance limit. First, while A^(n) interacts with terms 
which, in the large distance limit, become ijj K (x), the A\{n) interact with their 
derivatives d^ip K (x). We leave only the lowest order interaction terms, omitting 
the interaction terms containing the derivatives, therefore, the terms containing 
A^(x). Then, instead of a summation over all possible neighbours £ 6 Z 3 , we can 
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restrict the summation to the eight "non-trivial" basis lattice shifts uj £ {0,1} 3 : 
Terms which differ only by even lattice shifts become almost identical. This gives 

id t ip K (n) =(a l )i \ -i{i) K ,{n+ hi) - i/> K > {n - hi)) 
(35) +A-{n)^)i; Mn + hl)+ 2 ^'' {n - h ^ 

+ A%(n){r u )i^ K ,(n). 

Remember now that the operators f w are lattice approximations of the set of op- 
erators {1, 7 5 , 2Ii, 2Iij 5 } on the staggered lattice (see theorem [l]), and that these 
operators generate the Lie algebra of U(2)l x U(2)r — the left-right-symmetric ex- 



tension of the weak gauge group SU(2)l. Then, compare (35 1 with the continuous 
Dirac equation which interacts with some gauge field A^(x) via some representation 
f (1> e{l > 7 8 ,2/ i ,2I i 7 B }: 

(36) id t iP K {x) = (aTn(-^ + A?(x)t u )iI>(x) + A%{x)t^{x). 



We see that ( 35 ) is simply a discretization of ( 36 1 , which correctly takes into account 



that we have a staggered lattice. Remembering the lecture of fermion doubling, we 



should avoid premature claims that the continuous limit of pi) is (36 1. Nonetheless, 
for a given lattice equation (31 1, we can compute low energy effective fields A"(x) 



by averaging over the A^(n). And these effective fields A 1 ^ (x) interact with fermions 
like gauge fields with the gauge group U(2)l x U(2)r. This can be summarized in 
the following theorem: 



Theorem 2. A lattice deformation described by {3D defines effective low energy 
fields At(x). These fields interact with fermions in the same way as the gauge 
field of the left-right-symmetric extension U(2)l x U(2)r of the weak gauge group 
SU(2) L . 

Note here that it does not follow that for each continuous field configuration 
At (x) exists a physically meaningful deformed lattice giving this field in the con- 
tinuous limit. Instead, we obtain, below, two additional restrictions for the A^(x). 
In principle, we have not even proven here that there exist physically meaningful 
deformed lattices which give nontrivial fields A^(x) at all. Despite this, the con- 
struction of the A^(x) is not useless, but shifts the burden of argumentation: It 
is, now, the thesis that the fields A't (x) do not give the full group U(2)l X U(2)r, 
which requires justification. 

Note also that we do not claim here the existence of some effective gauge sym- 
metry. What we have found, is only a description in terms of effective fields A^ (x), 
which interact with the fermions in the same way as gauge fields. It is only the 
continuous limit of the interaction part which shows some gauge invariance. Differ- 
ent gauge-equivalent effective fields A^ (x) describe different geometric coefficients 
A|(n), thus, different deformed lattices. That means, on the fundamental level we 
have no gauge symmetry. Moreover, we have no theory about the equations of 
motion or a Lagrange formalism for these effective fields, nor on the fundamental 
level, nor in the continuous limit. This has to be left to future research. 

This theorem also allows the existence of other effective fields, which interact 
differently with the fermions. Here we have in mind the continuous limit of the 
fields A\(n), which give something like fields A^(x) interacting with the derivatives 
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of ip. Moreover, at the current state of research we cannot exclude, as well, that 
something similar to fermion doubling happens too, giving some other effective 



fields in the continuous limit. Nonetheless, the lattice equation (31 1 is already 
sufficient to restrict the interaction between these possible other fields and the 
fermions: Whatever these fields, the interaction is restricted to the maximal group 
U{2) L x U(2) R : 



Theorem 3. For lattice deformation described by the lattice equation 1 31 ), the 
resulting low-energy effective fields can interact with the fermions only via the set 
of operators f u G {1, 7 5 , 2ij, 2Iij 5 }, which generate the Lie algebra representation 
of the group £7(2) ^ x U(2)r. 

Indeed, the lattice coefficients interact with the fermion field ip only 

via the operators T£. But these operators become, in the continuous limit, 
t u G {1, 7 5 , 2/j7 5 }. And interaction with these generators cannot give more 
than an action of the group U(2)l x U(2)r. 

Now, it is reasonable to assume that the postulates [l][2]hold for deformed lattices 
as well as for undeformed ones. Thus, we require these postulates to our gauge- 
like correction terms as well. This leads to additional restrictions. While the 
preservation of the symplectic structure and rotational symmetry does not give 
anything new (11(2)^ x U(2)r is already unitary, preserves generations, and acts on 
all generations in the same way, thus, already preserves the symplectic structure as 
well as rotational symmetry), translational symmetry requires a further reduction of 
the maximal possible gauge group. Indeed, there has to be a preserved translational 
direction c in the leptonic sector. But the maximal group U(2)l x U(2)r does not 
leave any direction invariant. Thus, independent of the choice of the direction c, it 
cannot commute with translations. Any group, which commutes with translations, 
should be a nontrivial subgroup of U(2)l x U(2)r, and leave at least one particular 
direction c invariant. 

There are several possibilities for the translational direction c: ft may be left- 
handed, right-handed, or none of the above. The last case gives a more rigorous 
restriction of the group: The left-handed part 1 ~ 7 c as well as the right-handed 

part 1 ~t, 7 c would have to be preserved separately. Thus, to compute the maximal 
possible gauge group, we can ignore the last case. 

If the translational direction cis right-handed, we have the group U(2)l preserv- 
ing this direction, while U(2)r has to be reduced to some subgroup C/(l)i? C U (2)r. 
For this group, the charge of the translational direction c should be 0. Without 
restriction of generality, this charge can be taken as Ir = "t, 7 (J3 — 5), so that 
we obtain U(2)l x U(1)r as the maximal possible gauge group. Similarly, for 
left-handed c, we obtain the equivalent maximal gauge group U{2)r x U{1)l. 

Theorem 4. The maximal gauge group, which may be obtained from correction 
terms for lattice deformations, and is compatible with E(i) symmetry and symplec- 
tic structure, is a subgroup U(2) x U(l) of U(2)l x U{2)r. Without restriction 
of generality, this maximal group may be identified with U{2)l x U(1)r, gener- 
ated by the left-handed chiral U(2)l, and the right-handed chiral U(1)r with charge 

4.3. The EM field. At a first look, the EM field does not fit into any of the two 
classes of gauge-like lattice fields: ft acts nontrivially inside doublets, thus, is not 



16 



I. SCHMELZER 



a Wilson field. On the other hand, it has different charges on leptons and quarks, 
thus, does not fit into the our scheme for correction terms for lattice deformations. 

Nonetheless, we have already obtained it: It appears as a linear combination of 
the two types of gauge-like fields we have obtained. For the EM field we have 

(37) U{l) em C S{U(l) B x U(2) L x U(1) R ). 
Indeed, the EM charge 

1 1-t 5 1 1 + t 5 1 

(38) Q = 21 B + (h - i) = 2I B + -^(h i) + A^ih ~ g). 

is a combination of different charges: 1b comes from the diagonal of J7(3) c of Wilson 
gauge fields, J (I 3 — |) is part of U(2) l of weak interaction, and J (I3 — |) is 
the charge of U(1)r. A similar decomposition exists for the hypercharge: 

1 - T 5 1 + T 5 1 

(39) Y = 2I B - + A^ih i). 

4.4. The main result. Thus, all gauge fields of the SM are part of the maximal 
possible gauge group which can be constructed with our two ways to put gauge 
fields on the lattice: 

(40) C/(3) c x U(2) L x U(1) R D G S m = SU(3) C x SU(2) L x U(l) Y . 

Once these two types of lattice gauge fields are sufficient, let's add, in agreement 
with Ockham's razor, another postulate: 

Postulate 3. All gauge-like fields can be obtained as effective fields from the fol- 
lowing two types of lattice fields: 

• Wilson lattice gauge fields; 

• Correction terms for lattice deformations. 

With this postulate, the maximal possible gauge group is (without restriction 
of generality) U(3) c x U(2)l x U(1)r. It contains only three fields which are not 
part of the SM: The diagonal of U(3) c , the diagonal of U(2)l, and U(1)r. But, 
according to ( [39] ) , a linear combination of all three gives the hypercharge group 
U(1)y- Thus, we have only two additional gauge fields in the maximal possible 
gauge group. 

A further reduction can be obtained with the following postulate: 

Postulate 4. The gauge group should be a special group. 

This postulate can be justified using general properties of effective forces in con- 
densed matter theory. Volovik |27j writes: "An equilibrium homogeneous ground 
state of condensed matter has zero charge density, if charges interact via long range 
forces. For example, electroneutrality is the necessary property of bulk metals and 
superconductors; otherwise the vacuum energy of the system diverges faster than 
its volume. . . . The same argument can be applied to "Planck condensed matter" , 
and seems to work. The density of the electric charge of the Dirac sea is zero due 
to exact cancellation of electric charges of electrons, q e , and quarks, q u and qy, in 
the fermionic vacuum: Q vac — X^B<o(9 e + + = ®" Thus, to obtain neu- 
trality of the ground state, the trace of the charges has to be zero. Thus, the group 



In a similar way, a condensed matter approach promises to solve the cosmological constant 
problem 1271 : "The 3 He analogy suggests that a zero value of the cosmological term in the equilib- 
rium vacuum is dictated by the Planckian or trans-Planckian degrees of freedom: c9S V ac/fg M „ = 
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has to be a special group, with determinant 1. This postulate allows to get rid of 
the field U(1)b- This seems necessary, because already a very coarse consideration 
suggests that an U(1)b gauge field should be at least strongly suppressed: The 
Earth would be heavily charged, leading to a repulsive force on all baryonic matter. 
Thus, it should be weaker than gravity. A combination of gravity with this force 
would, moreover, lead to a violation of the EEP, thus, can be distinguished from 
pure gravity by tests of the EEP. Thus, an U(1)b field undetected by existing EEP 
tests should be much weaker than gravity, and, therefore, extremely suppressed in 
comparison with the EM field. 

Thus, there remains only a single additional gauge field. It can be chosen as an 
"upper axial gauge field" U(l)u, which acts on all upper particles (upper quarks 
and leptons) axially, with charge 

P 



(41) 



There is also a hyper- variant U(l)ft, which commutes with all SM gauge fields. It's 
charge is 



(42) 



i, 



u 



1 - 7 5 



I 



1 



R 



■r 



1+7 5 



2 " " ' 4 
This can be summarized in the following 



r 
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Theorem 5. The maximal possible gauge group G max with the following properties: 

• Postulate 1 ' preservation of the symplectic structure; 

• Postulate 2' preservation of E(3) symmetry; 

• Postulate^ realization on the lattice as 

— Wilson gauge fields; 

— Correction terms for lattice deformations; 

• Postulate restriction to special group; 

is S(U(3) C x U(2) L x U(l) R ) £ G SM x U(l) d . 

It acts, without restriction of generality, on the fermions with the standard action 
of Gsm and the U(l)g charge 



(43) 



1 + 7" 



h- 



2 



Now, the action of G ma x is anomalous. Thus, Gsm may be characterized as a 
maximal anomaly-free subgroup of G max - To obtain the SM, we could, as well, add 
anomaly freedom as an additional postulate: 

Postulate 5. The action of the group has to be anomaly-free. 

Because U (l)s defines also an anomalous extension of Gsm, w e can, in this case, 
even omit postulate |4j and obtain the following theorem: 

Theorem 6. The SM gauge group Gsm, with it's standard action, is a maximal 
group with the following properties: 



• Postulate 

• Postulate 

• Postulate 



preservation of the symplectic structure; 
2' preservation of E(3) symmetry; 
3' realization on the lattice as 



is the thermodynamic equilibrium condition for the "Planck condensed matter". Thus the equi- 
librium homogeneous vacuum does not gravitate. Deviations of the vacuum from its equilibrium 
can gravitate." 
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— Wilson gauge fields; 

— Correction terms for lattice deformations; 
• Postulate [3J anomaly freedom. 

But note that the standard rejection of gauge groups with anomalies is based 
on the canonical BRST approach to gauge field quantization, which requires exact 
gauge invariance. Instead, in our approach, gauge degrees of freedom are physical, 
and the Hilbert space is definite from the start, as discussed below in section [6j 
Therefore, the standard argumentation against anomalous gauge fields fails: With- 
out an indefinite Hilbert space structure, and without factorization of the (now 
physical) gauge degrees of freedom, there is no way to obtain a non-unitary theory. 
Therefore it seems premature to reject the extension U(l)u only because of the 
resulting gauge anomaly. So the author prefers, at the current state, theorem [5] 
and hopes for future observation of a new gauge field U(l)u- 

On the other hand, it seems also premature to claim that our approach predicts 
the additional field U(l)u- Indeed, we have to expect, that the anomaly leads to 
nontrivial effects during renormalization. As a result, the field U(l)u may as well 
become effectively suppressed. 

5. Fermion quantization 

Our approach to fermion quantization is a drastic departure from the common 
wisdom of quantum field theory. According to the standard approach, fermions do 
not have a classical configuration in the usual sense. Following Berezin [5J, instead 
of functions on a classical configuration space, we have, as some sort of classical 
limit, only Grassmann- valued fields. 

This is, obviously, incompatible with our geometric interpretation of fermion 
doublets as (C (£> A)(M 3 ), nor with the lattice model C(Z 3 ), which are classical, 
commuting, real fields, with a standard symplectic structure on the phase space. 
The appropriate way to quantize them would be canonical quantization. This 
seems, at a first look, impossible — last not least, we obtain in such a way only 
commuting operators. 

Despite this, we present here a way to obtain anticommuting fermion fields via 
canonical quantization. It consists of two parts, with a canonically quantized Z 2 - 
valucd field (spin field) as the intermediate step. The first part — to obtain a 
Z2-valued field from an M-valued field — is rather straightforward: All we need is a 
Z2-degenerated potential V(tp). The lowest energy states, then, define a Z2-valued 
field theory. 

The non-trivial step is from spin fields to fermion fields, or from commuting 
to anticommuting operators. Surprisingly, the lattice operator algebras appear to 
be isomorph. The isomorphism is, essentially, known in Clifford algebra theory. 
But this isomorphism is, first, nonlocal, and, second, not natural, but depends on 
some ordering between different lattice nodes. Moreover, it requires for dim > 1 
a nontrivial approximation of the Hamilton operator. Despite these caveats, the 
isomorphism between lattice fermion operators and Z 2 -valued "spin field" operators 
proves that fermions may be obtained via canonical quantization of real fields. 

In this paper, we need only canonical quantization of lattice theories with con- 
figuration space Q = M(Z 3 ) rcsp. Q = Z2(Z 3 ). But our considerations here do 
not depend on the dimension d = 3, so we consider here the more general case 
Q = M(Z d ) resp. Q — Z 2 (Z d ). Canonical quantization consists of the definition of 
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operators on the Hilbert space Jjf 2 (Q), and a Schrodinger equation 

(44) Wt*(q,t)=H*(q,t), q £ Q = ^(Z 3 ,Y), t£R. 

Thus, we always have continuous time. Note that in our condensed matter in- 
terpretation the lattice Z 3 is not a "discretization of space" ]R 3 itself. Instead, it 
enumerates elementary cells located in a continuous M 3 , where the state of the cells 
is described by some affine transformation Aff(3) C R 12 of M 3 . Nonetheless, in this 
section, there is no difference between our lattice Z 3 and a "discretization of space" 
Z 3 cl 3 , and our geometric interpretation of R 12 in terms of Aff(3) plays no role. 

5.1. Prom spin fields to fermion fields. Spin fields have the configuration space 
Z2(Z d ). On each lattice node n £ Z d we have the Pauli matrices a\ as operators: 

(45) a l n ai = 5 t] +ie l]k a k n . 
Spin field operators on different nodes commute: 

(46) [a^al] =2i5 mn £ijkV n - 

Instead, following Berezin[B], the fermion field operators - 0„, - 0* are usually con- 
sidered to be of completely different nature. They do not fit into the canonical 
scheme. Especially there is no configuration space Q. Operators related to differ- 
ent nodes do not commute. Instead, they anticommute: 

(47) = L< {CO = W>m,V>n} = 0. 

This difference seems to forbid an identification of fermions with spin fields. Despite 
this, the two operator algebras appear to be isomorph: 



Theorem 7. The operator algebra of spin field operators (45), (46 I and the alge 



bra of fermion operators (47 1 are isomorph. The isomorphism is not natural, but 
depends on the choice of an order > between the lattice nodes. 

To prove this, let's at first transform the operator algebras in each node into an 
equivalent form, by defining operators tp^: 

(48) ipl = i\) n + 1% = -*Wv, " <), 4% = 
This gives 

(49) ipi l ijji = Si j +iei jk ^, 



similar to (451. Now, for a given order > between the nodes of the lattice, the 
isomorphism is defined by the following formulas: 

(50a) V£ /a = II i>l = *l 

m>n 

(50b) <# a =v£ /a n^' 

Let's note that this isomorphism is well-known in the theory of Clifford algebras 
and allows to establish the isomorphism 

(51) Cl N ' N (R) Si AhiCf- 1 ^- 1 ^)) Si M 2 n(R). 

Indeed, the operators ipn an d ivpn generate (for a finite lattice with N nodes) the 
Clifford algebra Cl N ' N (M). On the other hand, M 2 n(R) is the operator algebra on 
the 2 Ar -dimensional space of Z 2 -valued functions on the same lattice. 
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Note also that (different from the a l n ) the operators ^\ do not act as local op- 
erators on the lattice. Instead, they act like on other nodes m > n. This is a 
necessary property of such an isomorphism, because, obviously, any local combina- 
tion of the commuting local operators a l n can give only to another set of commuting 
local operators. 

As a consequence, a Hamilton operator which "looks local" in terms of the i\) % n 
may appear nonlocal in terms of the a % n (which we consider to be "truly local" 
operators) and reverse. Fortunately, there are important examples of operators 
where this does not happen. First, there is the operator 

(52) #0 =~l^2 a n= \ Y. VnC 

n n 

Let's consider now operators with interactions between neighbour nodes. We 
are (for reasons which become obvious later) especially interested in the following 
linear combination: 

where hi are the d basic lattice shifts in the d-dimensional lattice 7L d . 

In the one-dimensional case, we have a natural (up to the sign) order >. For 
this order, we obtain: 

(54) h$ = i\ Y^Wnti+i + iM+i) = E ^n+i - CC+l 

n n 

Note that our operator itself is symmetric for spatial inversion n — > —n. But 
the representation in the asymmetric (in terms of the a l n ) operators ip^ hides this 
symmetry. 

5.2. The case of higher dimensions. Unfortunately, the transformation of the 
Hamilton operator in higher dimensions is not that simple. What we can obtain is 
only an approximation 

(55) « Hf =J2< + h i (^n +hi -r n r n+hi ) 



where 

{1 if n < n + hi 
-1 else 

The accuracy of this approximation obviously depends on the order >. Indeed, the 
error 



(57) ^^»=^ n G 



3 

n<m<n' 



resp. for cr^a^, , depends on the number and location of the nodes m located 
"between" (according to the chosen ordering) the "neighbour" (according to the 
lattice Z d ) nodes n, n' . Now, instead of the simple lexicographic order (which 
gives a™ +h . = 1) we propose to use another, more sophisticated order we name 
"alternating lexicographic order" . 

It has to be acknowledged that this order has been designed to give the result 
below. Fortunately, we can justify this choice of an order in another way: It gives a 
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Figure 2. The alternating lexicographic order 



better approximation of the original Hamiltonian operator, in the sense, that some 
algebraic properties of the original terms may be preserved exactly. 

Note that our interaction terms can be represented as a function of the differences 
of the operator a\ and its shift: 

(58) °\°\+ hi = 1-\{{1-T i )a 1 n )\ 

where is the shift operator on the lattice. This follows from (o^) 2 = 1 and the 
commutation relation [cr* , Tier*] = 0. Now, we propose to use an order which allows 
to preserve these properties exactly. That means, we want to replace the a\ by 
some a\ with exactly the same properties: 

(59) (a 1 n ) 2 = l,[^ n ,r l a 1 n } = 0, 
so that 

(60) aiai +hi « ^ +/t< = 1 - - t,)c£) 2 . 

For the simple lexicographic order, we have no way to define such a\. But it is 
possible for the alternating lexicogrpahic order. Let's define this order by induction. 
Let >k be the order defined for a k-dimensional lattice Z fe , and itk the projection on 
this lattice defined by the first k coordinates. Then we define >k+i by the following 
properties: 

• if n k+1 sg m k+1 then n sg fc+ i m; 

• else if rife+i(= rrik+i) is even and 7Tfcn ^ 7Tfcm then n ^k+i m; 

• else if rife + i(= rrik+i) is odd and iTkn ""few then n ^fe+i m. 

Thus, we use the inverse order inside the odd planes. Now the interaction term can 
be splitted in the following way: 

/C1\ 1/2 1/2 TT 3 -1/2-1/2 

n<im<n-\-hi 
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with 
(62) 



,V2 



n 



-1/2 



1/2 



n < m 
m, = n, 



n 

m < n + hi 
rtii — Hi + 1 



3 



and we obtain the properties ( 59 1 . The key is that for each node to with n < 
to, TOj = rii the shifted point to' = TjTO fulfils m' < n + h i7 thus, for each a^ n in the 
first term we find a corresponding cr^, in the second term. 

For our choice of >, the coefficients a™, fulfill the following relations: 



(63) 



n-\-2hi n n+hi -i 

a m+2hi> a n+hi a n+2hi — X ) 



5.3. Transformation of the lattice Dirac operator into staggered form. 

Now, the operator H — Ho + mHo appears to be a lattice Dirac operator. Indeed, 
let's consider the evolution equation defined by H: 



(64) 



(65) 



[H,il> n 



As a consequence of the relations ( 63 ) , the evolution equations (64 1, (65 1 give 



(66) 



where A 2 h is the lattice Laplace operator with doubled distance 2hi — a Laplace 
operator on a coarse lattice. 

The lattice Laplace operator A 2 h acts independently on 2 d different sublattices. 
Let's distinguish these sublattices by introduction of 2 d different lattice functions 
enumerated by elements of k = (aci, . . . , k^) G {0, Using the denotation *ip n = 
tp*, we define 



(67) 



j,Kl + ... + Kd 



V', 



on 



= k mod 2. 



Each of the 2 d lattice functions ip K (n) is defined on a "coarse lattice" containing 
the nodes of type Uj = 2n^ + and lattice spacing 2hi. Now, the lattice operator 
A 2 h acts as the simple Laplace operator on each of the 2 d functions ip K (n). In the 
continuous limit, each ip K (n) gives a function ip K (x) which fulfills the Klein-Gordon 
equation 



(68) 



d 2 Mx, t) = 0? - ™ 2 )M%, t) = o. 



The lattice Dirac equations (64 1, (65 1 now establish a connection between these 



2 d lattice fields. We can define now 2 d x 2 d matrices (a 1 )* , so that the original 



lattice equations (64), (65 1 transform into 



id t ip K (n) = [H,ip K (n)] 
(69) = ^-jfa^'^ln + fti) - ^(n-fei)) +m^> K -(n) 



on n = n mod 2. Note that, because of the factor * k i+---+ k <j j n (67l, equation (69) 
connects only the fields ip K (n), and it's adjoint only the ^p* K {n). 
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This equation is our lattice Dirac equation (18 1 on the staggered lattice (20 1, 
but already in its quantized form, with anticommuting fermion operators ip K (n), 
and with a mass term. 

5.4. Prom spin fields to scalar fields. Spin fields are already a much more 
classical object in comparison with the original fermion fields. But for our approach 
we need even more classical objects, namely real-valued fields. 

Fortunately, this is not problematic at all. We can embed the spin field as 
an effective description of the lowest energy states of a scalar field with a Z 2 - 
symmetric potential which gives two Z 2 -symmetric vacuum states. For example, 
we can consider t/? 4 theory in M. d with negative mass parameter /1 2 : 

(70) Sf=^{B t <p) 2 -{d i <p) 9 )-V{<p)withV{<p) = "' ' X 




The two minima of the potential are ip(x) — ±tpo with ipo = y^v- 

If the system is near (po, it is convenient to use the c-variable <j(x) = ip(x) — ipo 
so that 

(71) V(a)^ 1 2 (2^W + y [ j^ + ^a\ 

This describes a scalar field with mass v2/^ and some interactions. 

Instead, we are interested only in the lowest energy states of this theory. Let's 
consider at first the simple case of dimension d — 0, where QFT reduces to ordinary 
quantum theory. If we have energies much below /i, only the two vacuum states 
1 &±(< r ) with (^il^l^i) ~ ±<po are important. But the true eigenstates of energy 
are 

(72) tf 0/1 (p) = i(*+(p)± 

Between them, we have an energy gap of order 

r<pu ^3 



rvo „ u 
(73) A = Ex - E ~ exp(- / y/V(<p) - E d<p) ~ exp(-^). 

Jo 



With increasing /i the mass of the a field increases, but the energy gap A decreases 
exponentially Without any conspiracy, this leads to two different domains: a high 
energy domain, with energies of order /i, and a low energy domain, with energies 
of order A, where the whole theory reduces to the two-dimensional space spanned 
by ^0/1 • Reduction to this subspace gives 

(74a) (p — > fo^ 1 with <p = I v]/ . ip^^ip 



(74b) 7r = — -» 7r er 2 with tt q = I ^ • d^xdip 



1 , j TT E + Ei 



(74c) H -> H - -Acr J with H = , ( 

For dimension d > 0, at least as long as the momentum k is sufficiently small, 
we have a similar situation for each of the modes <p(x) — exp(ikx)ip. For suffi- 
ciently large /1, and sufficntly low energies under consideration, the theory reduces 
to an effective theory where we have only two degrees of freedom for each mode. 
Effectively, the configuration space reduces from ^(Z^M) to J r (Z d ,Z 2 ). 
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Last not least, let's consider typical lattice theory interaction terms which may 



appear in the reduction for a Lagrangian of type ( 70 1 . We consider lattice ap- 
proximations where only neighbour nodes have nontrivial interaction terms. Let 
n,n' = n + hi be these neighbour nodes, d = 1. One possibility is to use \ (ir n + ir n / ) 
to approximate n(x) and ^(<^n — fn') to approximate diip(x). Then, the reduction 
gives an effective Hamiltonian 

(75) H = 1 -{tt 2 + c 2 (d x ^) 2 + V{ip)) -> c + c x a l n a l n , + c 2 a 2 n a 2 n , + c 3 a 3 n 

for some constants c,-. The lattice Dirac operator corresponds to ci = —c 2 = 1, 
thus, can be obtained in this scheme. 

As a consequence of this quantization method for fermions, we obtain some 
analogon of a "supersymmetric partner" of the fermions. This partner can be very 
heavy without any conspiracy. At the current state of research, no indications 
about their masses can be given. 

5.5. Conclusion. Given that our proposal to obtain fermionic fields via canonical 
quantization of real fields is a very drastic departure from conventional wisdom 
about fermion quantization, let's summarize here the results of this section: 

• The connection between real fields and spin fields is conventional and un- 
problematic. The only unconventional clement we have used here is to use 
a discretization in space but not time. 

• In one-dimensional space, for our choice of the order, the Hamilton operator 
does not require any approximation. Thus, for spatial dimension 1, the 
problem of obtaining fermions from canonically quantized Z2-field operators 
a l n on the lattice is solved exactly and completely. Once the question, if 
fermions may be obtained via canonical quantization of real fields, at least 
at a first look seems to be dimension-independent, this exact result for 
dimension 1 seems to be a powerful argument for a positive answer. 

• The isomorphism between the operator algebra of spin field operators a % n 
and algebra of fermion operators tp n , tpn i s > m every dimension, and for every 
choice of the order, an exact isomorphism. It is only the Hamilton operator, 
which has to be approximated in higher dimensions, and which depends on 
the choice of the order. But the definition of the quantum operator algebra 
is, certainly, the more important, fundamental part of quantization. Thus, 
the most important part of the quantization scheme is exact too. 

A better justification of the choice of the order, as well as a evaluation of the error 
made during the approximation, have to be left to future research. Nonetheless, 
in the light of the exact results we have obtained, the general problem, if fermion 
fields may be obtained via canonical quantization of real fields, may be considered 
as solved. 



6. Some notes about gauge field quantization 

After the fermion quantization, the approach to quantum gauge fields is a second 
drastic departure from conventional wisdom: With the correction coefficients for 
lattice deformations, we propose to use, essentially, a non-gauge-invariant regular- 
ization for the chiral gauge fields. Moreover, we reject, in general (that means, even 
for the Wilson lattice gauge fields) the concept of BRST quantization in its current 
form. 
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In our condensed matter approach, the gauge fields have to be obtained as ef- 
fective fields from more fundamental fields. In the case of the Wilson gauge fields, 
these will be degrees of freedom which describe the material between the cells, which 
effectively influences the connection between two neighbour cells in such a way, that 
it has to be compensated by the gauge transformation corresponding to the edge. 
Instead, the correction coefficients for lattice deformations are, essentially, degrees 
of freedom of the configuration of the cells. Thus, they have been already quantized 
by the quantization of the fermionic part. In above cases, the gauge fields have not 
to be quantized themself, in a separate quantization procedure, but to be derived, 
as effective quantum fields, similar to phonons, from the more fundamental quan- 
tum theory)^] How to do this has to be left to future research. Despite this, let's 
give here answers to some objections against our approach, which are based on the 
standard BRST approach to gauge field quantization. 

A first objection against our construction of weak gauge fields in section |4.2| is 
that it presents a lattice regularization for chiral gauge field theory. But to obtain 
such a regularization is a famous problem of chiral lattice gauge theory |15j . and 
there are various no-go theorems for such regularizations. But the regularization 
problem of chiral gauge theory is the problem to find a gauge-invariant regular- 
ization. Our regularization has no exact gauge invariance on the lattice. Instead, 
we have only approximate gauge invariance — the generators of the gauge group 
are associated with nontrivial lattice shifts. Thus, our regularization is not in con- 
tradiction with the various no-go theorems for regularizations with exact gauge 
invariance. 

This answer leads, in a natural way, to a second objection. Last not least, 
people have tried to find a regularization with exact gauge invariance not just for 
fun, but for a good reason — to quantize chiral gauge fields. The problem is that 
the standard procedure to quantize gauge fields — BRST quantization — depends 
essentially on exact gauge invariance of the theory. Without exact gauge invariance, 
it fails miserably: What remains is a non-unitary theory. But this failure is a special 
problem of the manifestly Lorentz-covariant Gupta-Blcuer approach to gauge field 
quantization, which starts with an indefinite Hilbert space structure. Following 
a proposal of Gupta [14] and Bleuer [7] for QED, manifest relativistic invariance 
is reached in the BRST approach using an unphysical "big space" with indefinite 
Hilbert metric. A physical interpretation of this big Hilbert space would lead to 
negative probabilities, which is nonsensical. To get rid of the states with negative 
probability, restriction to an invariant subspace and factorization is used. But these 
operations depend on exact gauge invariance. If gauge invariance fails, the result 
is fatal for the whole approach. 

But there is even a classical alternative — the earlier approach of Fermi [TTj 
and Dirac 10J, where the Hilbert space is definite, but Lorentz covariance is not 
explicit. Whatever may go wrong, the Hilbert space remains definite, and at least 
a probability interpretation of the results is possible. Another possibility is our 
approach, where quantization is done in terms of a more fundamental theory, which 

^ Let's note in this connection the approach of Volovik 1271 , who uses examples of usual 
condensed matter, for example 3 He — A, as analogies to fundamental physics. Especially, he 
considers the appearance of effective gauge fields and gravitational fields in usual materials. His 
observations appear close to our approach. For example |27| : "The 3 He — A analogy also suggests 
that all the degrees of freedom, bosonic and fermionic, can come from the initial (bare) fermionic 
degrees of freedom." 
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is based, of course, on a definite Hilbert space. Especially, in our model we have 
quantized the degrees of freedom of the cells (whose positions define the lattice, 
and, therefore, also the correction coefficients for lattice deformations) with classical 
canonical quantization. Then, the gauge fields are derived as effective fields, similar 
to phonons. In this approach, there is no place for an indefinite Hilbert space 
structure, and, therefore, a non-unitary theory, to appear. Let's note that, once 
we don't require exact gauge symmetry, the standard rejection of anomalous gauge 
fields is no longer justified. 

The rejection of the manifestly Lorentz-covariant approach leads to another ob- 
jection: It has not been introduced without reason too. Manifest Lorentz covari- 
ance, on one hand, simplifies computations. This is, obviously, not a decisive argu- 
ment. More serious is that, in our approach, we do not have relativistic invariance. 
Indeed, our construction from the start violates Lorentz covariance, and in many 
different ways: First, we handle time and space in different ways, having a lattice 
only in space. Then, even a spacetime lattice would violate the symmetry of the 
continuous limit. Moreover, the operators cr^, for spatial spinor rotations on our 
staggered lattice, are nonlocal, and, therefore, do not define an exact representa- 
tion of the algebra su(2) An approach, which violates Lorentz invariance on the 
fundamental level, has to explain, how it will be recovered in the large distance 
limit . 

Fortunately, this question has been, at least partially, addressed by the deriva- 
tion of the Einstein equivalence principle (which includes local Lorentz covariance) 
in our theory of gravity [25] . Because of the importance of this question, we give 
an introduction into this theory in appendix [XJ On the other hand, the connection 
between the postulates of this theory of gravity and the approach described here 
remains unclear. Especially, our fermion quantization procedure does not lead to 
a unique "speed of light" for different fermions. For the solution of this question, 
a result of Chadha and Nielsen [8] may become important: They considered mass- 
less electrodynamics with different metrics for the left-handed and right-handed 
fermions; their model thus violates Lorentz invariance. They found that the two 
metrics converge to a single one as the energy is lowered. Thus, in the low-energy 
corner Lorentz invariance becomes better and better. 

7. Discussion 

Many questions have to be left to future research. This includes: 

• Symmetry breaking; 

• The search for a Hamilton operator for a general configuration of cells, 
which would allow the derivation for other regular crystallographic lattices 
as well as for lattices with deformations and defects; 

• The large distance limit, especially renormalization group equations; 

• The connection between our lattice model for the SM and the theory of 
gravity. They are conceptually and metaphysically compatible, but math- 
ematically yet unrelated; In the large distance limit of our cellular lattice, 
we will obviously have notions like density p, average velocity v l and some 
stress tensor cr u , which follow continuity and Eulcr equations. Nonetheless, 



'This violation of Lorentz invariance appears also for four-dimensional staggered fermions. 
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the open question is how to obtain the corresponding Lagrange formalism, 
which is postulated in the theory of gravity (see (|82|). 

Especially symmetry breaking promises to be interesting. First, we need it. The 
E(3) action does not define a symmetry of the SM: The fermionic mass terms clearly 
violate the rotational symmetry between the three generations. Moreover, the whole 
construction of section [5] which creates effective Z^-fields from the original M-fields, 
violates translational symmetry: We cannot add real constants to Z2-valued fields. 
Thus, to obtain the Hamiltonian of the SM, even to obtain fermion fields at all, we 
have to break E(3) symmetry. 

On the other hand, some points suggest large differences between the SM sym- 
metry breaking scheme and the symmetry breaking we need in our approach: 

• There is an interesting correspondence between gauge invariance of the SM 
and gauge invariance in our lattice model: Strong interactions have exact 
gauge symmetry, in our model as well as in the SM, weak interactions 
not: Nor in our model, nor in the SM after symmetry breaking, where the 
masses of the gauge bosons as well as the masses of the fcrmions violate 
gauge symmetry. From point of view of this correspondence, it would be 
strange if, in a first step, we obtain exact gauge invariance for weak gauge 
fields, which, then, in a second step, becomes broken again. Thus, even 
if the symmetry breaking participates in giving gauge bosons mass, the 
starting point will, probably, not be a gauge-invariant theory. 

• A standard argument for a separate Higgs sector is that spatial isotropy 
is unbroken. But in our model, spatial isotropy is broken by the mass 
terms. Thus, this argument in favour of a separate Higgs sector fails. The 
question appears if we need a separate Higgs sector at all. The application 
of Ockham's razor — don't introduce Higgs fields without necessity - 
suggests that there will be no Higgs sector at all. 

Thus, we need symmetry breaking, but for very different reasons, and of a dif- 
ferent symmetry. Therefore we expect the symmetry breaking in our model to be 
very different from the SM scheme. The only common thing will be, possibly, the 
general idea of symmetry breaking, and that it has to give fermions mass. 

Despite these open questions, our cellular lattice model, being of surprising sim- 
plicity, already has a place for all standard model particles observed until now. 
The fermionic content of the SM is predicted exactly, with all three generations. 
Already this part, taken alone, gives our model large explanatory power: Even for 
small natural modifications (two or four instead of three generations, two or four 
instead of three colors, clcctroweak singlets, triplets or quartets instead of doublets) 
where would be no corresponding condensed matter model of similar simplicity. 

Then, the maximal gauge group compatible with few simple postulates gives 
the SM gauge group almost exactly, leaving room for only one additional gauge 
field U(l)if. This is also a quite impressive result, given that the maximal pos- 
sible unitary gauge group on the fourty eight Weyl fermions of the SM — the 
48 2 -dimensional group £/(48) — has been reduced to a 13-dimensional subgroup 
S(U(3) C x U(2)l x U(1)r), which already contains the 12-dimensional gauge group 
SU(3) C x SU{2) L x U(1) Y of the SM. 

The comparison of these results with those of other approaches, like string theory, 
will be left to the reader. 



28 



I. SCHMELZER 



Last not least, the model is conceptually compatible with a metric theory of 
gravity with GR limit. Thus, it gives also hope for a development of a "theory of 
everything" , which unifies the SM fermions and gauge fields with gravity. 



Appendix A. Gravity 

For metric theories of gravity there is a simple way to obtain a condensed matter 
interpretation, closely related to the ADM decomposition [T] or the geometrody- 
namic interpretation |28j . The preferred frame of such an interpretation defines an 
ADM decomposition of the four-metric g a p into a scalar field, a three-vector and a 
definite three-metric. We identify these fields with density p, velocity v l and stress 
tensor a 13 of some form of condensed matter in the following way: 

(76a) g ot) V=g = P, 

(76b) g 0i ^=pv i , 

(76c) g l3 V^9 = Prfv j - a ij . 

For these condensed matter fields, we would like to have continuity and Euler 
equations: 

(77a) d t p+d l (pv l ) = 

(77b) d t (pv l ) + 0<(/wV - a ij ) = 0. 



They coincide with the harmonic conditions for the metric ( 76 ) : 

(78) d a (g af, V=g) = 0. 

This condensed matter interpretation is, therefore, possible for all metric theories 
of gravity, which include the harmonic condition as a physical equation. A simple 
theory with this property is general relativity in harmonic gauge. One variant of 
GR in harmonic gauge is to add a non-covariant term to the GR Lagrangian which 
enforces the harmonic conditions: 

(79) L = Z a g aa v^g- + L GR (g a P,i J matter ) 

for some constants S Q . Its dependence on the preferred coordinates X a (x) can be 
made explicit Q 



(80) L = -^S^XlX^-g + L GR (g a P, ^ matter ) 

This explicit form is useful because it allows variation over the preferred co- 
ordinates. We have to take care — the four functions X a (x) + SX a (x) have to 
define a valid system of coordinates — but nonetheless variation is possible and 
gives Euler-Lagrange equations for the X a of the same form as for usual fields. We 
obtain: 



The dependence of some expression on the preferred coordinates X a (x) is, by definition, 
explicit, if, after a formal replacement of occurrence of X a (x) by four scalar fields Y a (x), the 
resulting expression is covariant. So, in a a replacement X a (x) — > Y a (x) changes nothing, 
the resulting expression a is not covariant, thus, the dependence on the preferred coordinates 
is implicit. Instead, in the form a^X®^, the replacement gives the expression a^Y^, which is 
covariant, because Y° is considered as a scalar field. 
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(81) 



SS 
SJCy 



a/3 



?a^ 7 ), 



thus, the preferred coordinates X a (x) are harmonic. The Lagrangian (79 1 obvi- 
ously defines a metric theory of gravity with Einstein equivalence principle. In the 
limit S Q — > we obtain the Einstein equations. The terms g aa \f—g do not depend 
on partial derivatives of the metric, therefore the limit S Q — > is natural for small 
distances and weak fields. 



But, as long as we simply postulate the Lagrangian (80), we have no explana- 



tion for these properties. The classical argument against the Lorentz ether may 
be raised: It needs some conspiracy, does not give an explanation for relativistic 
symmetry Is it possible to derive this Lagrangian from some postulates which are 
more natural for a theory with preferred frame? 



Theorem 8. The Lagrangian {80) follows from the following two conditions: 



(82a) 
(82b) 



5S - (» 



SX° 



) + dj(pv z v j 



of (81 ) is defined by a particular solution (given by the first, non-covariant term of 



Indeed, given ( 76 1 , the equations ( 82 I are equivalent to ( 81 1 . The general solution 



(80 1) and the general solution of the homogeneous problem 



(83) 



SS 
SX" 



0, 



thus, modulo a covariant Lagrangian. The covariance of the Lagrangian is what 
we take here as the definition of the Lagrangian Lqr of general relativitjj^] in ( 80 ) 

□ . 



Now, to postulate the equations (82) docs not require much conspiracy. Instead, 



they can be seen as a combination of the Noether theorem with the standard in- 
terpretation of continuity and Euler equations, as conservation laws for energy and 
momentum in condensed matter theories. Indeed, if the Lagrangian, in its explicit 
form, has a symmetry X a — ► X a + c, the Euler-Lagrange equation for the preferred 
coordinates does not depend on the X a themself, but only on its partial deriva- 
tives. In this case, the Euler-Lagrange equation automatically obtains the form of 
a conservation law: 



(84) 



SS 
SX" 



-d 



dL 



dX 







Thus, the left hand side of (82 1 define the Noether conservation laws related with 
translation in time and space. On the right hand side we have the continuity and 

^Notc that wo use here the most general understanding of general relativity, where the Einstein- 
Hilbert Lagrangian is only the lowest order term, and higher order terms, or terms with higher 
order derivatives in the metric, are, in principle, allowed, as long as they are covariant. This 
understanding is standard for effective field theories — in the large distance limit, only the lowest 
order terms survive. 
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Euler equations — the conservation laws for energy and momentum in condensed 
matter theory. To identify left and right hand sides is a very natural postulate for 
a condensed matter theory. 

More details and consequences of this theory of gravity can be found in |25j . 
Especially, the gauge-breaking term stops (for the correct sign of the constants) the 
black hole collapse and prevents the big bang singularity. The condensed matter 
approach to gravity solves many quantization problems of canonical GR quantiza- 
tion: The notorious "problem of time" [T7] simply disappears. Together with the 
black hole collapse the related information loss problem [2~3] disappears too. 

Appendix B. The Dirac operator on the exterior bundle 

We give here the basic formulas for the Dirac operator in the exterior bundle 
(see, for example, [22]) ■ The exterior bundle or de Rham complex A = J2k=o^ k 



^fc=o ■ 

consists of skew-symmetric tensor fields of type (0, fc), < k < d, which are usually 
written as differential forms 

(85) ip — ipi 1 ...i h dx il A • • • A dx ik 6 A fe 

The exterior bundle A has dimension 2 d in the d-dimensional space. The most 
important operation on A is the external derivative d : A k — > A fc+1 defined by 

where i q denotes that the index i q has been omitted. It's main property is d 2 = 0. 
In the presence of a metric, we have also the important ★-operator A k — > A d ~ k : 

1 

fc!" 

with -k 2 = (— l) k ( d ~ k >sgn(g) . This allows to define a global inner product by 



(87) (*V0h + x...<„ = ,., jr n ■ ■■g' lk3k V\g\^ 1 . 



(0,7/0 = J </>A(*ip) = J V> a (*</>) 

It turns out that the adjoint operator of d* : A k — > A k ~ 1 of d is 

(89) d* = (-l) kd+d+1 * d* 

Note that the expressions for * 2 and d* depend on the order of the form fc, which 
is not nice. But a minor redefinition of the 7k- operator allows to solve this problem. 
For the operator 

(90) * = i k(d - k h 
the resulting expressions no longer depend on fc: 

(91) * 2 = sgn(s), d* = (-i) d+1 * d* 

In this general context we can define the Laplace operator as 

(92) A = dd*+d*d. 

Then, the Dirac operator (as it's square root) can be defined as 

(93) D = d+d*, 

so that A = D 2 . Indeed, we have d 2 — as well as (d*) 2 = 0. 
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B.l. Discretization of the Dirac operator. The special geometric nature of 
the exterior bundle allows to define a nice doubler-free discretization of the Dirac 
equation on a general cell complex. Such a cell complex consists of cells Ci of 
dimension k = dim(ci), which are embeddings of the fc-dimensional unit cube I k 
into the manifold so that the image of the boundary is part of the image of lower- 
dimensional cells of the complex, and the image of all cells of the cell complex covers 
the whole manifold. 

On such a cell complex, fc-dimcnsional differential forms are represented on the 
lattice by their integrals over the fc-dimensional cells Cj of the cell complex: 



The external derivative defines in a similar natural way a derivative for functions 
on the mesh, with the same most important exact property d 2 = 0. 

For the definition of the ^-operator we need a metric and a dual mesh. A metric 
on the manifold defines in a natural way for every cell Ci it's area aj = a(ci) > 0. 
For a triangulation on Euclidean background, the values oii depend on each other. 
But in the general case they may be considered as independent variables, which 
approximate the metric on the cell complex. In the following we consider them as 
given and defining the metric. 

A dual mesh is a mesh with cells Cj so that for each cell Ci of the original 
mesh with dimension k we have a corresponding "dual cell" c; of dimension d — k 
which intersects only the cell Cj, in a single point, orthogonally and with positive 
intersection index. The metric defines the areas a, of the dual cells in a similar way. 
Now, the lattice Hodge ^-operator may be defined as 



and maps functions on the mesh to functions on the dual mesh. Note that the 
dual of the dual mesh has the same cells as the original mesh, but possibly with 
different orientations. Therefore, for * 2 we obtain an additional factor ( — \^ k ( d ~ k ) 
as in the continuous case. 

Thus, we can define the exterior derivative as well as the Hodge * operator on 
the lattice preserving their algebraic properties d 2 = 0, * 2 = (— l) fc ( d_fc ). As a 
consequence, the remaining part of the theory can be transferred on the lattice too. 

It is a general consequence of the geometric character of the continuous Dirac 
equation as well as its lattice discretization that the lattice discretization does not 
have doublers. See, for example, [21 HI 151 124"] . 
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